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Abstract 

SU(2) lattice gauge theory with four flavors of quarks is simulated at nonzero 

chemical potential \x and temperature T and the results are compared to the 

predictions of Effective Lagrangians. Simulations on 16 4 lattices indicate that 

at zero T the theory experiences a second order phase transition to a diquark 

condensate state. Several methods of analysis, including equation of state 

fits suggested by Chiral Perturbation Theory, suggest that mean-field scaling 

describes this critical point. Nonzero T and \x are studied on 12 3 x 6 lattices. 

For low T, increasing \x takes the system through a line of second order phase 

transitions to a diquark condensed phase. Increasing T at high fj,, the system 

passes through a line of first order transitions from the diquark phase to the 

quark-gluon plasma phase. Metastability is found in the vicinity of the first 

order line. There is a tricritical point along this line of transitions whose 

position is consistent with theoretical predictions. 
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I. INTRODUCTION 



Recently there has been a resurgence of interest in QCD at nonzero chemical potential 
for quark number. Arguments based on instantons and phenomenological gap equations 
support the old idea that diquark condensation and a color superconductivity phase 
transition [[|^| will occur at a critical chemical potential somewhat greater than that at 
which nuclear matter forms. Nuclear matter is expected to be the ground state of QCD above 
a chemical potential hn, which is slightly less than one-third of the proton's mass. These 
arguments neglect the forces of confinement and may only be reliable at asymptotically large 
chemical potentials. Since diquarks carry color in real QCD, a quantitative understanding 
of confinement and screening is needed to estimate the energies of the states which control 
the phases of the system. Unfortunately, brute force lattice simulation methods, which are 
so useful at T ^ and fi = 0, do not yet provide a reliable simulation algorithm for these 
environments in which the fermion determinant becomes complex. 

Given these problems, theorists have turned to simpler models which can be analyzed 
and simulated by current methods. One of the more interesting is the color SU(2) version of 
QCD which addresses some of the issues of interest Jl|], 0, ||. In this model diquarks do not 
carry color, so their condensation does not break color symmetry dynamically. Instead, the 
diquark condensed phase resembles a superfluid rather than a superconductor. The critical 
chemical potential is also expected to be one-half the mass of the lightest meson, the pion, 
because quarks and anti-quarks reside in equivalent representations of the SU (2) color group. 
Hence hadron flavor multiplets include both mesons and diquarks. Chiral Lagrangians can 
be used to study the diquark condensation transition in this model because the critical 
chemical potential vanishes in the chiral limit, and the model has a Goldstone realization of 
the spontaneously broken quark- number symmetry |6|-|9| . The problem has also been studied 
within a Random Matrix Model at non-zero /i and T [[K]]. Lattice simulations of the model 
are also possible because the fermion determinant is real and non-negative for all chemical 
potentials. One hopes that these developments will uncover generic phenomena that will also 
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apply to QCD at nonzero chemical potential. Furthermore, most of the theories that aim at 
a description of the phase diagram of 3-color QCD at nonzero baryon chemical potential and 
temperature can be implemented to describe the phase diagram of 2-color QCD at nonzero 
baryon chemical potential and temperature. Therefore, an important motivation to conduct 
Lattice simulations is to check the validity of these theories in the 2-color case. If a theory is 
not correct in the two-color case, it is very doubtful that the predictions of the corresponding 
3-color theory can be trusted. 

Preliminary lattice simulations of the SU(2) model with four species of quarks, simulation 
data and an Effective Lagrangian analysis of aspects of the T-jj, phase diagram were recently 
published |TT| , [12 l. It is the purpose of this article to present the continuation of that work 



on larger lattices, closer to the theory's continuum limit with a focus on the system's phase 
diagram at nonzero /i and T. Very early work on this model at finite T and /i was performed 
by |13[ . A simulation study of the spectroscopy of the light bosonic modes will be presented 



elsewhere (n| . 

In our exploratory study |T2[ of this model's phase diagram, we found an unanticipated 
result. At fixed but small quark mass, which insures us that the pion has a nonzero mass 
and chiral symmetry is explicitly broken, we found that at high T and \i there is a line 
of transitions separating a quark-gluon plasma phase from one with a diquark condensate. 
Along this line there is a tricritical point, labeled 2 in Fig. 1, where the transition switches 
from being second order (and well described by mean field theory) at relatively low /i to a first 
order transition at an intermediate /i value. We will present evidence here for metastability 
along the line of high fi and T. We have argued elsewhere ]12| that this simulation result, 



the existence of a tricritical point 2, has a natural explanation in the context of chiral 
Lagrangians. Following the formalism of || we argued that trilinear couplings among the 
low lying boson fields of the Lagrangian become more significant as /i and T increase and 
they can cause the transition to become first order at a /i value in the vicinity of the results 
found in the simulation. Using hindsight, this behavior should not have come as a surprise. 
fj, plays the role of a second 'temperature' in this theory in that it is a parameter that 



controls diquark condensation, but does not explicitly break the (quark-number) symmetry. 
The existence of two competing 'temperatures' is a characteristic of systems which exhibit 
tricritical behavior. A tricritical point is also found in Chiral Perturbation Theory, as well 



as in a new Random Matrix Theory model [|15|,|nj. 

Schematic Phase Diagram of Diquark Condensation in the Temperature-Chemical Potential Plane 
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FIG. 1. Schematic Phase Diagram in the T-jj, Plane. The thin(thick) line consists of sec- 
ond(first) order transitions. The dashed line denotes a crossover. Point I labels a critical point 
and Point 2 labels a tricritical point. The existence and position of point 1 is a matter of conjecture 
in the two color model. 

The conjectured phase diagram of Fig. [I] also has a critical point labeled 1 which connects 
the line of first order transitions to the dashed line of crossover phenomena extending to the 
/i = axis where we expect a finite T chiral transition between conventional hadronic matter 
and the quark-gluon plasma. For sufficiently light quarks in the four flavor SU(2) model, 
the /j, = transition is known to be first order both from theoretical arguments [[T^IOl and 



simulations |19j. For quark mass m = 0.05 we find that the transition is smoothed out to a 
crossover region as fi is increased until the dynamics favors diquark condensation. Clearly, 
the results will depend on the precise value of the quark mass m. At m = 0.05 we find 
that the transition at fi = 0.10 is as smooth as that at fi = 0.0, while the transition at 
\i = 0.20 is noticeably sharper. This result suggests that if m were chosen somewhat smaller 
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than m = 0.05, the scenario with the critical point 1 and the first order line might emerge. 
However, this critical point 1 might be absent so the crossover line would intersect the curve 
separating the diquark condensation phase from the normal phase. We will be pursuing this 
suggestion in separate simulations. The simulations done here are not accurate enough to 
probe the details of the region of the phase diagram Fig. |T] where the two first order lines 
split apart. Simulations at lower quark masses on larger lattices might be needed. 



In QCD at nonzero baryon chemical potential Fodor and Katz |2(| have estimated the 



position of the critical point 1 using a modified, multivariable Glasgow algorithm pi] . It 
will be interesting to find point 1 in the four flavor SU(2) model at smaller m values using 
our conventional simulation methods and then checking whether the Fodor-Katz method 
can find it with competitive accuracy. If point 1 is absent, the Fodor-Katz method should 
be able to find the intersection of the crossover line with the curve delimiting the diquark 
phase. 

We also present extensive and accurate simulation results on a 16 4 lattice, essentially 
at vanishing T. The critical indices of the transition at nonzero fi are measured and favor 
mean-field scaling in agreement with Chiral Perturbation Theory, analyzed through one loop 
corrections, and simulations of 2-color QCD in the strong coupling limit [fSJ. The critical 



indices f3 mag and 5 are emphasized. Equation of State fits, suggested by Chiral Perturbation 
Theory, also support these conclusions. As a third method to determine the critical indices, 
we use the static scaling hypothesis. This analysis is also compatible with mean field, 
although this approach is unable to predict the index 5 accurately because of its sensitivity 
to estimates of the critical chemical potential, /i c . Unfortunately, because simulations at 
small A are extremely costly in computing requirements, we have not been able to approach 
the critical point closely enough to make definitive measurements, insensitive to the detailed 
nature of the critical scaling. In a recent paper in which we study quenched versions of 
this and closely related theories, we have indicated how sensitive measurements of critical 
scaling are to the form of the scaling variable |]23[| . This paper does indicate, however, that 
there are classes of scaling variables for which there exist extended regions outside of the 



true scaling window, where the order parameter scales with the correct critical exponents. 



II. SIMULATION METHOD AND THE ALGORITHM 

The lattice action of the staggered fermion version of this theory is: 

Sf = £ \x[P(jfi + m) X + 1\[X T T2X + XT2X T ]\ (1) 

sites ^ A ' 

where the chemical potential \i is introduced by multiplying links in the +t direction by 
and those in the —t direction by e _At The diquark source term (Majorana mass term) 
is added to allow us to observe spontaneous breakdown of quark-number on a finite lattice. 
The parameter A and the usual mass term m control the amount of explicit symmetry 
breaking in the lattice action. We will be particularly interested in small values of A and 
the extrapolation A — > for a given m to produce an interesting, realistic physical situation. 
This is the system that has been studied analytically using effective Lagrangians at non-zero 
chemical potential \i H0], and non-zero temperature |L5|. 
Integrating out the fermion fields in Eq.l gives: 



pfaffian 



Ar 2 A 
-A T \t 2 



det(^U + A 2 ) (2) 



where 

A = p(fj) + m (3) 
Note that the pfaffian is strictly positive, so that we can use the hybrid molecular dy- 



namics 



25| method to simulate this theory using "noisy" fermions to take the square root, 
giving Nf = 4. 

For A = 0, m / 0, /i ^ we expect no spontaneous symmetry breaking for small \i. 
For /i large enough (fi > according to most approaches including Chiral Perturbation 

Theory H-fl) we expect spontaneous breakdown of quark number and one Goldstone boson 



- a scalar diquark. (The reader should consult || for a full discussion of the symmetries 
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of the lattice action, remarks about spectroscopy, Goldstone as well as pseudo-Goldstone 
bosons, and for early simulations of the 8 flavor theory at A = 0.) 



III. SIMULATION RESULTS AND ANALYSIS 



We now consider the simulation results for the Nf = 4 theory on 16 4 , 8 4 and 12 3 x 6 
lattices, measuring the chiral ((xx)) an d diquark condensates ((x T T2x)), the fermion number 
density, the Wilson/Polyakov line, etc. 
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FIG. 2. (x Tt 2X) as a function of \x on an 8 4 lattice at (3 = 1.5, m = 0.025. The unlabeled 
squares are a linear extrapolation to A = 
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Let us first briefly mention our 8 4 simulations. These were performed at (3 = 1.5, a 
relatively strong coupling, and m = 0.025. These extend our earlier work with the same 
P but m = 0.1. Using a smaller quark mass increases the validity of chiral perturbation 
analyses. In addition, we are performing simulations on a 12 3 x 24 lattice with the same 
parameters to examine the spectrum of Goldstone and pseudo-Goldstone bosons of this 
theory. For this purpose it is necessary to have the pion mass significantly lighter than that 
of non-Goldstone, non-pseudo- Goldstone excitations. We have simulated at A = 0.0025, and 
0.005. Our measurements of the diquark condensate are given in figure ^. Each 'data' point 
represents a 2000 molecular-dynamics time-unit simulation. The fact that a simple linear 
extrapolation of our measurements to A = is very small, not only at /x = 0, where we know 
that the condensate vanishes, but also at /x = 0.1 leads us to believe that the condensate 
vanishes in this limit, probably up to /x ~ 0.2. (Note that is is known that there is some 



small but finite value of /x, below which this condensate vanishes p2| .) For /x much greater 
than 0.2 but less than the saturation value it is reasonably clear that the condensate does 
not vanish in the limit as A — > (on an infinite lattice). Hence there is a phase transition 
for some value of /x to a state where quark number is spontaneously broken by a diquark 
condensate. The transition occurs over a relatively small range of /x, but does not appear 
to be first order. Because of the relatively small scaling window this implies, we will not be 
able to analyze critical scaling until we perform measurements at more /x values within this 
scaling window. 

The chiral condensate remains approximately constant up to /x = /x c after which it falls 
rapidly approaching zero for large /x. The fermion density rises slowly and approximately 
linearly from zero for /x > /x c . Above /x ~ 0.6, it increases more rapidly approaching its 
saturation value of 2 at /x ~ 1.0. 

We now turn to measurements on a 16 4 , 'zero temperature' lattice. We simulated the 
SU(2) model at a relatively weak coupling (3 = 1.85, within the theory's apparent scaling 
window, in order to attempt to make contact with the theory's continuum limit. The quark 



mass was m = 0.05, as in ||12| , and a series of simulations were performed at A = 0.0025, 



0.005, and 0.01 so that our results could be extrapolated to vanishing diquark source, A = 0. 
Simple linear and spline extrapolations were used and these procedures appeared to be 
sensible away from the transition. None of the conclusions to be drawn here will depend 
strongly on the limit. In fact, everything we say could be gathered from our data at a 
fixed (small) A value, 0.005 or 0.0025, say. However, only in the limit of vanishing A do we 
expect real diquark phase transitions (at least when the transitions are second order), so 
it is particularly interesting and relevant to investigate the A — > limit. The quantitative 
agreement between fits at fixed, small A values and those using extrapolated data enforce 
the prevalent idea that simple extrapolation procedures are reliable except in the immediate 
vicinity of the critical point where non-analyticities are inevitable. 

Because we are forced by practical considerations to work at A values large enough that 
the condensate is relatively smooth over the transition region, we have to deduce critical 
scaling from the curvature of the condensate further from fi c than we would like. Thus 
we run the risk that we are outside the scaling window and the apparent scaling we are 
seeing does not represent true critical scaling. However, we saw for quenched theories [^J at 
relatively weak couplings, that the form of the natural scaling variable is such that it has an 
inflection point when expressed in terms of the simplest scaling variable /x — fi c . This leads 
to an extended domain where the order parameter scales with the true critical exponents. 
It is reasonable to assume that dynamical theories behave similarly. 

In Table I we present the raw data on the 16 4 lattice at the three A values. Table II 
exhibits the results of extrapolating this data to A = 0. Typically, one thousand time units 
of the Hybrid Molecular Dynamics algorithm were accumulated to generate these data sets. 
The error bars in the data sets of Table I account for the correlations in the raw data sets. 
The low A runs were our most CPU intensive. 

Since the Hybrid algorithm suffers from systematic errors proportional to the square 
of the discrete time step used in integrating its stochastic differential equations forward 
in Monte Carlo time PS[, we were forced to a small discrete time step. The simulations 



reported here used time steps as small as dt = 0.0025 to control these errors. Even then, 
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because our statistical errors are so small, these systematic errors might not be negligible. 
We will later indicate cases where we used dt values which were larger than optimal and 
have seen signs of such errors. 

Now let us discuss some results. In Fig. |3] we show the diquark condensate linearly 
extrapolated to A = plotted against the chemical potential /i. The quark mass m was 
fixed at 0.05, the coupling was (3 = 1.85 and the linear extrapolation used the 
A = 0.0025 and 0.005. 

Diquark Condensate on 16 A 4 lattice, beta=1.85, m=0.05, lambda=0.00(linear) 
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FIG. 3. Diquark Condensate vs. /x 

We see evidence to a quark-number violating second order phase transition in this figure. 
The dashed line is a power law fit from jj, = 0.30 to 0.35 which predicts the critical chemical 
potential of // c = 0.2860(2). The power law fit is good, its confidence level is 48 percent, and 
its critical index is (3 mag = 0.54(3) which is consistent with the mean field result (5 ma g = 1/2, 
predicted by chiral perturbation theory , J?j including one loop corrections. Note that the 
quark mass is fixed at m — 0.05 throughout this simulation, so chiral symmetry is explicitly 
broken and the transition of Fig. 2 is due to quark number breaking alone. 

In figure § we show the diquark data extrapolated to A = quadratically using the raw 
data at A = 0.010, 0.005 and 0.0025. The fit to the data at /i ranging from 0.30 through 
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0.35 is also shown. The fit predicts a critical point \x c = 0.2870(2), has a confidence level 
of 47 percent and has a critical exponent of I3 mag = 0.54(4), again in agreement with mean 
field theory. 

Diquark Condensate on 16 A 4 lattice, beta=1.85, m=0.05, lambda=0.00(spline) 
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FIG. 4. Diquark Condensate vs. /x 
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Further evidence for j3 ma g = 1/2 can be obtained directly from the raw data without 
any extrapolations. For example, in figure |5] we show the A = 0.0025 data and a power law 
fit ranging form fj, = 0.30 through 0.38. We find fi c = 0.2749(2), (3 mag = 0.54(5) with a 
confidence level of 24 percent. 

Diquark Condensate on 16 A 4 lattice, beta=1.85, m=0.05, lambda=0.0025 
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FIG. 5. Diquark Condensate vs. p. 

We also considered extrapolations of the form (x Tt 2X) = P where p is the solution of an 
equation of the form p 3 = a p + bX with a and b determined as functions of p by performing 
an independent fit at each p. This form is suggested my mean-field scaling forms, which 
would predict that a should be proportional to p — p c close to the transition. However, these 
fits were so poor that this approach was abandoned. 

In addition, we considered the fermion number density which is predicted to rise linearly 
from p c and found good agreement with this expectation, as shown in figure |6|. 

Fermion Number Density on 16 A 4 lattice, beta=1.85, m=0.05, lambda=0.0025 
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FIG. 6. Fermion Number Density vs. p. 

The fit predicts the critical point p c = 0.2765(5), with a critical index 1.17(7) and a 
confidence level of 13.1 percent. Plots and fits of the fermion number density at the other 
A values were equally good. In fact the fermion number density shows relatively little A 
dependence. This is in agreement with Chiral Perturbation Theory 0. 

Finally, we considered the A = 0.005 simulations and extended the measurements up to 
large p — > 1.2, shown in figure 0. We see the curious phenomenon observed in the past: as 
p becomes large, the condensate falls to zero. The point is that once p becomes of order 
unity, the discreteness of the lattice comes into play and suppresses the density of states 
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which forces the condensate to fall toward zero. The quark-number density approaches 
2 (per lattice site), the maximum allowed by Fermi statistics. This is a lattice artifact. 
Nonetheless, we can try a power law fit to the diquark condensate over the entire region 
of \x from 0.30 to 0.60 where the curve is increasing. The result is shown in the figure and 
the fit predicts f3 mag = 0.49(7) and \i c = 0.2651(2). Since x 2 /dof ~ 5, the fit is poor, but 
the deviations from this fit are small enough that this fit should be considered a reasonable 
zeroth order approximation. 

Diquark Condensate on 16 A 4 lattice, beta=1.85, m=0.05, lambda=0.005 
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FIG. 7. Diquark Condensate vs. \i. 



In summary, j3 mag measurements are relatively easy and self-consistent. For all the 
parameters and procedures we tried, we found good agreement with mean field theory. It 
was important, however, to be working at relatively weak coupling, (3 = 1.85, in the theory's 
scaling window where such an extended scaling window exists. Experience has shown that 
at stronger coupling, the true scaling is described by effective Lagrangians closer to the 
non-linear sigma model form, where the scaling window in \i — \i c is small and the extended 
scaling regime exhibits scaling in /i — // c but with an apparent critical exponent (3 m which is 
half the true value |2g^J2§. 

A second critical index in which we are interested is 5, which controls the response of 
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the system to symmetry breaking fields at the critical point /i c , 

(x T r 2X ) = AX 1/S (4) 

in the limit of small A. Mean field theory predicts 5 = 3. This scaling law Eq. 4 is usually 
difficult to use in practice because first it requires accurate estimates of /x c , and second it 
generally suffers from a small scaling window in A. Not having a precise estimate of fi c is 
probably our biggest problem, since as we have shown in |23||, fits to the extended scaling 



regime can underestimate /i c by ~ 8%. In order to obtain some useful data, however, we 
simulated the 16 4 lattice for four estimates of fi c each over a range of A from 0.002 to 0.010. 
First we note that the chosen dt for the A = 0.002 simulations was too large and although 
plotted on our graphs, these points were excluded from our fits. 

We considered 0.2855, 0.2870, 0.2920 and 0.2970 as estimates of fi e and tried fits of the 
form A\ l / & + B. The data and fits are shown in the next four figures, which we discuss in 
turn. 

The /i = 0.2855 data predicted 5 = 4.0(6) with a confidence level of 91 percent, but with 
a large, negative constant B = —0.17(7). The fit considers the data from A ranging from 
0.004 to 0.010 because the data at A = 0.002 has large dt 2 errors. These results are shown 
in figure |] The fact that B is large and negative suggests that this estimate for /i c is too 
low. 

In figure |9] we consider the same ideas for the estimate \i c = 0.2870 and now find 5 = 
2.5(2) and B = —0.02(2), but with a rather poor confidence level of a few percent. For /i c = 



0.2920 (figure [Hi]), we find 5 = 2.1(2), but B is distinctly different from zero, B = 0.05(2). 



The confidence level is, however, very good, 96 percent. Finally, for // c = 0.2970 we find 



5 = 2.7(6) and B = 0.02(4) with a confidence level of about one percent. See figure [TT. 
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Diquark Condensate on 16 A 4 lattice, chemical potential = 0.2855, m = 0.05 
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FIG. 8. Diquark Condensate vs. A for ji c = 0.2855 



Diquark Condensate on 16 A 4 lattice, chemical potential = 0.2870, m = 0.05 
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FIG. 9. Diquark Condensate vs. A for ji c = 0.2870 
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Diquark Condensate on 16 A 4 lattice, chemical potential = 0.2920, m = 0.05 
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FIG. 10. Diquark Condensate vs. A for fi c = 0.2920 



Diquark Condensate on 16 A 4 lattice, chemical potential = 0.2970, m = 0.05 
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FIG. 11. Diquark Condensate vs. A for fi c = 0.2970 



In summary, 5 is in the mean field 'ballpark', but the uncertainly in the simulation results 
is large. The results give us no reason to doubt the application of mean field theory because 
they are all in the vicinity of 5 = 3 with some fits higher and some lower. The sensitivity of 
the simulation estimates to fi c is considerable. We learn from these measurements that only 
a focused set of measurements that determine \i c very well will produce a useful estimate of 
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8. 



IV. COMPARISON WITH CHIRAL PERTURBATION THEORY 

It is interesting to compare the lattice data to the theoretical predictions of Chiral Per- 
turbation Theory. Two-color QCD at non-zero baryon/quark- number chemical potential 
has been thoroughly studied at the tree level and at next-to-leading order within Chiral 
Perturbation Theory [|6|-|8|. At next-to-leading order, it was found that the phase transi- 
tion between the normal phase and the superfluid diquark phase is second order, that the 
critical chemical potential is given by fi c = and that the critical exponents are given 

by mean-field theory. Furthermore, based on the structure of the loop corrections, it was 
conjectured that the critical exponents are given by mean field at any (finite) order in Chiral 
Perturbation Theory ||. 

Overall, the next-to-leading-order corrections do not dramatically change the leading- 
order results ||. The corrections turn out to be small, as long as \x is close to /i c . Therefore 
it should be enough to compare the lattice data to Chiral Perturbation Theory at leading- 
order, at least for the lattices sizes, couplings and parameters considered here. 

In Chiral Perturbation Theory at leading order, the diquark condensate at non-zero 
diquark source is given by 



where (xx)o is the chiral condensate at zero chemical potential, zero temperature, and zero 
diquark source, and a is given implicitly by 



In the equation above tan</> — X/m, and m, is the mass of the Goldstone modes at chemical 
potential fi = 0, temperature T = 0, and diquark source A = [0. Furthermore the 
quark- ant iquark (chiral) condensate is given by 



(X T ^2X) = (xx)osina, 



(5) 



4/i 2 cos a sin a = m 2 sin(a — <p) . 



(6) 
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(xx) = (xx)ocosa. (7) 

At nonzero diquark source there is a crossover between the normal phase and the diquark 
condensation phase. However, Chiral Perturbation Theory at leading order predicts that 
(Ax) 2 + (x Tt 2X) 2 does not depend on the chemical potential, the diquark source, or the quark 
mass. This relation does not hold at next-to- leading order in Chiral Perturbation Theory || . 
However, for small chemical potential, quark mass and diquark source, Chiral Perturbation 
Theory at leading order should give an accurate description of these observables. Therefore, 
we expect that (xx) 2 + (x Tt 2X) 2 depends slightly on the quark mass and the diquark source, 
in the same way as the quark-antiquark condensate depends on the quark mass in three-color 
QCD at zero temperature and chemical potential [p7fl . The behavior of (xx) 2 + (x Tr 2X) 2 is 



presented in Figure [12]. These graphs confirm what was expected from Chiral Perturbation 
Theory: (Xx) 2 + (x Tr 2X) 2 is close to constant for small enough fi and A. Notice that it is 
impossible to see the critical chemical potential, which is around 0.29 as determined in the 
previous section, by looking at (Xx) 2 + (x Tr 2X) 2 - The points at A = 0.010 show the effects 
of finite dt 2 errors. The points which lie higher were generated in simulations with dt = 0.01 
while the 4 lower points were generated in simulations with dt = 0.005. 



(xx) 2 + (x t t 2 x) 

0.3- 



. 2 



o.i- 



0.25 0.3 0.35 0.4 
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FIG. 12. (xx) 2 + (x Tt 2X) 2 f° r different diquark sources and chemical potentials. The pluses 
are the lattice data for A = 0.025, the xs for A = 0.004, the triangles for A = 0.005, the squares for 
A = 0.006, the diamonds for A = 0.008, and the dots for A = 0.010. a) shows that (xx) 2 + (x Tr 2X) 2 
is almost constant for 0.28 < \i < 0.30 and for a given diquark source, b) shows the transition 
region on an expanded scale. 

We could compare the diquark condensate obtained on the lattice with the leading-order 
result (H). But we notice that within Chiral Perturbation Theory at leading order we have 
that 

\xx) 

This observable can also be used as an order parameter of the diquark condensation phase. 
In Chiral Perturbation Theory at leading order, for a given quark mass, diquark source and 
chemical potential, the ratio @ only depends on the value of the pion mass m n through the 
minimization equation (j^). It is therefore a very suitable observable to use in the comparison 
between the lattice results and Chiral Perturbation Theory. 

Hence we fit the lattice result for (x Tr 2X)/ (xx) with Chiral Perturbation Theory 
This is a one-parameter fit. Since we use Chiral Perturbation Theory at leading order to 
analyze our results, we can only take into account data corresponding to chemical poten- 
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tials and diquark sources that are small enough so that (xx) 2 + (x Tr 2X) 2 is constant for a 
given diquark source. The previous figure shows that we have to limit our analysis to data 
corresponding to 0.28 < \i < 0.30, where (xx) 2 + (x Tt 2X) 2 is close to constant for a given 
diquark source. If we use the 11 points that correspond to A = 0.004, 0.005, and 0.006, and 
0.28 < // < 0.30, we get that fi c = 0.3027(1). This is an acceptable fit since x 2 /dof = 2.3 
(remembering that the norm of the condensate is only approximately constant over this 



range.) This one-parameter fit and the data we used to perform it are shown in Figure [13|a. 
If we now use the 18 points that correspond to A = 0.0025, 0.004, 0.005, 0.006, and 0.008, 
and 0.28 < ft < 0.30, we get that /i c = 0.3028(1), and the fit has a x 2 /dof = 3.1. If we 
use the 25 points that correspond to A = 0.0025, 0.004, 0.005, 0.006, 0.008, and 0.010, and 
0.28 < fi < 0.30, we get that [i c = 0.3026(2). The fit is worse than the previous ones with 
a x 2 /dof = 12. This one-parameter fit and the data we used to perform it are shown in 
Figure p~3]a. Again we note that the 3 points at A = 0.010 which lie above the curve are 
those with dt = 0.01 which have larger dt 2 errors, and are largely responsible for this large 
X 2 - 

Noting that the new estimates for [i c are significantly higher than our previous estimates, 
we consider the inclusion of the 'data' at /i = 0.31 and = 0.32 in our fits. If we include the 
measurements from = 0.31 to the 25 points that correspond to A = 0.0025, 0.004, 0.005, 
0.006, 0.008 used in the fits of the previous paragraph, the x 2 /dof increases from 3.1 to 4.9. 
Including both the /i = 0.31 and [i = 0.32 measurements further increases the x 2 /dof to 5.9. 
In both cases the new /i c is entirely consistent with the estimates of the previous paragraph. 
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FIG. 13. Comparison between Chiral Perturbation Theory and the lattice results for 
(x T T 2X) / (xx) ■ The pluses are the lattice data for A = 0.025, the xs for A = 0.004, the trian- 
gles for A = 0.005, the squares for A = 0.006, the diamonds for A = 0.008, and the dots are 
for A = 0.010. The results of two different fits with Chiral Perturbation Theory and the data 
we used for these fits are depicted in the solid curves for the different diquark sources: a) for 
A = 0.004, 0.005, 0.006 only; b) for all As. 
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In summary, we find that Chiral Perturbation Theory at leading order describes the data 
well, with a critical chemical potential given by fi c = 0.3027(3). This is somewhat higher 
than our earlier estimates in the previous section, and might explain why our 5 estimates 
were less than compelling. It is also understandable in terms of what was learned from the 
quenched studies performed in In Chiral Perturbation Theory, the critical exponents 
are given by mean-field theory even at next-to-leading order. From the quality of the fits of 
the lattice results with Chiral Perturbation Theory that are presented above, we conclude 
that the critical exponents measured on the lattice are consistent with mean field. 

We notice that Chiral Perturbation Theory at leading order does not describe the data 
very well at large chemical potential and diquark source. This is presumably a sign that 
higher order corrections have to be taken into account. The introduction of /x breaks the 
original symmetry and allows the condensates to vary independently. This can, in part, be 
implemented by replacing the effective potential based on non-linear sigma models, by a phe- 
nomenological effective potential suggested by linear sigma models, where the condensates 
are less tightly constrained [E3 . 



V. SCALING FUNCTION 

The static scaling hypothesis states that the data for (x Tt 2X) /(/■* — n c ) l3mag for different 
diquark sources and chemical potentials should collapse onto a single curve for /j > \i c when 
plotted against A/(/x — /i c ) A [p8fl - The gap exponent A = f3 mag 5; A = 3/2 in mean-field 
theory. We can fit the lattice results to the scaling function if we assume some form for the 
scaling function. We will try several possibilities for the form of the scaling function: linear, 
of the form f(x) = a + bx s , of the form f(x) = a + bx + cx 2 ^ 3 and the simplest mean- field 
form. The number of parameters changes according to the assumption for the form of the 
scaling function. We present our results in Table 4 using part of the data (/i = 0.297, 0.30, 
0.31, 0.32, 0.33, and 0.34). 

The different forms assumed for the scaling function give similar results for the critical 
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chemical potential and for the two critical exponents that we studied. The best fit corre- 
sponds to a critical chemical potential that is consistent with the result we found in the 
previous section, /i c = 0.2948(7), a value of (3 mag = 0.58(4) which is consistent with mean 
field, and a value of 5 = 2.28(3) which is about three quarters of the mean-field result. None 
of these fits are very good: the best fit has a x 2 /dof = 17. It corresponds to the scaling 
function of the form f(x) = a + bx + cx 2 / 3 , as presented in Table VI. In Figure [14|, we show 
the scaling function together with the 95%-confidence ellipse for the parameters (3 mag and S 
using the best fit. 
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0.25 0.5 0.75 1 

mag 

FIG. 14. Scaling function and 95%-confidence ellipse for the critical exponents [3 mag and 8 
from our best fit, which corresponds to fi c = 0.2948, (3 mag = 0.58, and 5 = 2.28. The pluses are 
the lattice data for A = 0.025, the xs for A = 0.004, the triangles for A = 0.005, the squares for 
A = 0.006, the diamonds for A = 0.008, and the dots for A = 0.010. 

The quality of these fits is rather poor. Furthermore, we find that there is a sizeable 
correlation between the different fit parameters: \i c is strongly anti-correlated with f3 mag and 
mildly with 5, whereas j3 mag and 5 are mildly correlated with each other. Therefore, the 
results from the fits using the scaling function are somewhat suspect. This can be readily 
seen by using the same data as above plus the four points that correspond to /i — 0.2920. 
The best fit we get has a x 2 /dof = 66. It gives fi c = 0.269(8), (3 mag = 0.7(1), and 5 = 1.4(2). 
Notice that /i c is much lower than before, that the increase in f3 mag is sizeable, and that 5 
has changed by almost 40%. This result illustrates the problems usually encountered in the 
use of the scaling function. Since the scaling function is in general different above and below 
the critical point, the choice of data is crucial, and it is difficult to get a clear result from 
this type of analysis. The comparison with Chiral Perturbation Theory presented in the 
previous section does not suffer from these problems. 

Finally if we impose mean-field exponents, using the same part of the data as above (// = 



24 



0.297, 0.30, 0.31, 0.32, 0.33, and 0.34), the best fit we get has poor quality: x 2 /dof = 243. 
It gives fi c = 0.297(8), which is compatible with the results obtained in the previous section 
where we compared the data with Chiral Perturbation Theory. If we try to use scaling 
functions derived from the mean- field equation of state as in P3"| , the quality of the fit is 
marginally better, with a y 2 /dof = 205, and it gives /i c = 0.308(2). 



VI. RUNS AT FINITE TEMPERATURE 

We now present our 12 3 x 6 simulations to investigate the interior of the phase diagram 
of Fig. m . We are particularly interested in scanning the diagram starting at fi = and 
following phase transitions and lines of crossover phenomena at high T out to large fi. 
In particular, the low mass four flavor theory should have a first order transition on the 
/i = axis and this transition should penetrate into the phase diagram. In addition, the 
T = transition to a diquark condensate at \i c should penetrate into the phase diagram and 
separate the normal phase from the phase with a diquark condensate as shown in Fig. 1. It 
would be interesting to understand if and how these two lines of transitions merge inside the 
phase diagram. Finally, we know from past simulations and Effective Lagrangian analyses 
that there is a line of first order transitions at high /i and high T which separates the 



diquark condensate phase from the quark-gluon plasma phase ||12|| . We want to confirm the 
first order character of the high /i part of this line and to find the low /i end of this line 
where it changes to second order at a tricritical point, labeled 2 in Fig. |l|. This point has 
also been found within Chiral Perturbation at nonzero temperature and chemical potential 



0. 



We are hopeful that some of these features of the phase diagram will also occur in other 
models, so it would be useful to confirm them in this relatively simple setting. Perhaps, 
QCD at nonzero Baryon Number chemical potential has a.T-fi phase diagram with some of 
these features. 

We begin by considering the fi = axis. Our simulations are at nonzero fermion 
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mass m = 0.05 and non-zero diquark source A = 0.005 (equivalent to running at m = 
V0.05 2 + 0.005 2 ). m = 0.05 is large enough that we expect the transition between hadronic 
matter and the quark-gluon plasma to be weakened and perhaps softened into a crossover. 
This is, in fact, what we find. (At lower quark mass we would expect to find a first order 
transition.) In Fig. |1| we show data for the chiral condensate and the Wilson Line and 
confirm a smooth but rapid crossover in the vicinity of (3 = 1.9. 

Wilson Line and Condensates on 12x12x12x6 lattice, chem pot=0.00, m=.05, lambda=0.005 
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FIG. 15. Wilson Line and Condensates vs. (3 for [i = 0.00. 



Measurements at fi = 0.10 gave similar conclusions: there is a clear but smooth crossover 



in the vicinity of /3 — 1.9, as shown in figure [16 
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Wilson Line and Condensates on 12x12x12x6 lattice, chem pot=0.10, m=.05, lambda=0.005 
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FIG. 16. Wilson Line and Condensates vs. [3 for [i = 0.10. 



It is interesting that a plot of the same quantities at fi = 0.20 displays a noticeably 
sharper crossover. Figure |T7] suggests that fi = 0.20 is near a critical point, such as the 
point 1 in our generic phase diagram, Fig. |I|. At this rather large quark mass m = 0.05, this 
appears merely to be a more pronounced crossover but, perhaps, if m were chosen smaller 
than 0.05, an actual line of first order transitions would be found for the region corresponding 
to \i > 0.20 in this phase diagram. 

Wilson Line and Condensates on 12x12x12x6 lattice, chem pot=0.20, m=.05, lambda=0.005 
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FIG. 17. Wilson Line and Condensates vs. (3 for fi = 0.20. 
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Once we increase fi to 0.25, the induced diquark condensate becomes more significant at 



low T, as shown in figure 18 



Wilson Line, Condensates on 12x12x12x6 lattice, chem pot=0.25, m=.05, lambda=0.005 




FIG. 18. Wilson Line, and Condensates vs. [3. 



At \i = 0.30 the diquark condensate becomes the most interesting and rapidly varying 
quantity. The chiral condensate is relatively smooth and changes from 0.493(6) at (3 = 1.7 to 
0.217(1) at (3 = 1.9 while the diquark condensate has varied from 0.352(5) to 0.0400(3) over 
the same range of coupling. The Wilson Line varies from 0.057(3) to 0.241(2) over the same 
range. As shown in the figure [T9], the diquark condensate experiences a sharp transition 
between (3 = 1.80 and 1.90. 



28 



Wilson Line and Condensates on 12x12x12x6 lattice, chem pot=0.30, m=.05, lambda=0.005 
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FIG. 19. Wilson Line and Condensates vs. (3. 



That transition becomes even sharper as we increase fx to 0.35 as shown in figure 20 

Wilson Line and Condensates on 12x12x12x6 lattice, chem pot=0.35, m=.05, lambda=0.005 
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FIG. 20. Wilson Line and Condensates vs. [3. 
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At n = 0.40, the simulations indicate a first order transition near (3 = 1.85(3). In figure ^1 
we show both the diquark condensates and the Wilson Line and see strong suggestions of 
discontinuities. 
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Wilson Line and Condensates on 12x12x12x6 lattice, chem pot=0.40, m=.05, lambda=0.005 
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FIG. 21. Wilson Line and Condensates vs. (3. 

The really solid evidence for a first order transition comes from the time evolution of the 
observables at (3 = 1.87 which show signs of metastability. For example, in the figure we 
show the time evolution of the diquark condensate and display several tunnelings between 
a state having a condensate near 0.15 and another with a condensate near 0.40. 

Evolution of Diquark Condensate at beta=1.87, chemical potential=0.40, m=.05, lambda=0.005 
0.7 
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FIG. 22. Diquark Condensate vs. Computer Time. 
Runs at = 0.50 and 0.60 show clear discontinuities in both the diquark condensate and 
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the Wilson Line near (3 = 1.87. In figure |23| we show the results for \x = 0.60. 



Wilson Line and Condensates on 12x12x12x6 lattice, chem pot=0.60, m=.05, lambda=0.005 
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FIG. 23. Wilson Line and Condensates vs. j3. 
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Finally, we scanned the phase diagram at low T and variable fi. We accumulated data 
at — 1.30, relatively strong coupling and near vanishing T, on the 12 3 x 6 lattice and 
took measurements over a wide range of /i, from 0.15 to 0.60. The results for the diquark 
condensate are shown in the figure |24] which shows a transition near fi = 0.25. The power 
law fit shown in the figure produces a critical index (3 mag = 0.29(4) at a critical chemical 
potential /i c = 0.2453(2). The quality of the fit, which extends from [i = 0.26 to 0.38, is 
rather poor, having a confidence level of only a fraction of one percent. As we have seen in 
quenched simulations and in simulations of QCD at finite chemical potential for isospin such 
small estimates of the critical index are often an indication that the scaling is best described 
in terms of the scaling form for the effective Lagrangian in the non-linear sigma-model class. 
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Diquark Condensate on 12x12x12x6 lattice, beta=1.3, m=.05, lambda=0.005 
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FIG. 24. Diquark Condensate vs. /i. 



The estimate for the critical chemical potential, fi c = 0.2453(2), meshes well with our 
measurements at high T above, which indicated that the diquark condensate only started 
to be numerically significant for fi > 0.25. 

The continuous phase transition is also seen in the induced fermion-number density 
shown in the figure |25|. The dotted line fit to that data, which extends all the way from 
fi = 0.26 to 0.56, has a confidence level of 42 percent, a critical index 0.99(2) and an 
estimate for the critical chemical potential, \i c = .2282(2). The scaling law obtained here 
is in good agreement with field theory expectations of unity Note that such linearity for 
these relatively strong couplings has a wider range of validity than would be indicated by 
the scaling window, as is seen in the explicit scaling forms from effective Lagrangians ||[7|]- 
Since we only have a single A value, we have not attempted a fit to such a form. 
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Fermion Density on 12x12x12x6 lattice, beta=1.3, m=.05, lambda=0.005 
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FIG. 25. Fermion-Number Density vs. fi. 



VII. CONCLUSION AND OUTLOOK 

In this article we have presented a study of the phase diagram of 2-color QCD at nonzero 
baryon/quark-number chemical potential and nonzero temperature. We have found two 
phases: a "normal" phase where quark number is conserved and a phase with a diquark 
condensate which spontaneously breaks quark number. At fi = 0, for smaller quark masses, 
there should be a first order finite temperature "deconfinement" transition which divides 
the hadronic from the quark-gluon plasma phase. Since, even in this case we expect the 
line of first order transitions emerging from this point to terminate at small fi, beyond 
which there should be a rapid crossover but no transition, the hadronic and quark-gluon 
plasma phases are not distinct. This is the phase we call the "normal" phase. The chiral 
condensate is non-zero everywhere. However, where it is large the system is more hadronic 
in nature; where it is small, the system is more plasma-like. In the "diquark" phase, the 
chiral condensate decreases with increasing diquark condensate, approaching zero for large 
\i (even before saturation sets in). We have found that the phase transition between the 
"normal" phase and the diquark phase for small T is second order and compatible with mean 
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field theory. These simulations therefore agree with the predictions of Chiral Perturbation 
Theory through next-to- leading order ||. 

We have also found that the second order phase transition between the "normal" phase 
and the diquark phase becomes first order when \i and T are increased. Therefore the \x- 
T phase diagram contains a tricritical point. It corresponds to a chemical potential near 
two-thirds of the pion mass, and a temperature only slightly lower than the transition 
temperature between the hadronic phase and quark-gluon plasma phase at fi = 0. The 
tricritical point is also found in Chiral Perturbation Theory at nonzero temperature and 
chemical potential [15|]. Therefore, this phase diagram is consistent with the one obtained 



within Chiral Perturbation Theory [ |i~5|j . 

Now that the phase diagram of this theory has been mapped out, we can turn to more 
quantitative properties. In particular, we plan on analyzing the gauge configurations for their 
topological content. There are interesting predictions [[HJ based on Effective Lagrangians 
for the size, density and interactions among the instantons of the model at large \x which 
can be investigated by cooling methods which have been very successful at vanishing [i. In 
addition, we can simulate the phase diagram at low /i and high T with lighter quarks and 
search for the critical point 1. The SU(3) analog of this point is thought to be accessible to 
heavy ion collisions planned for BNL's RHIC. 

Furthermore, the spectroscopy of the model, with an emphasis on Goldstone and pseudo- 
Goldstone bosons, is also under consideration. Since the theory's light modes control its 
critical behavior and thermodynamics, quantitative results on spectroscopy are quite impor- 
tant. 

There are also interesting effective Lagrangian predictions for QCD with chemical po- 
tentials associated with the light quarks |3T|-|33|j . We are studying QCD with a chemical 



potential associated with isospin f2q , |T5|| . This model's phases at nonzero fi and temperature 
T are expected to be very similar to those discussed here ||31|| . Although we cannot attack 
the SU(3) theory with a large Baryon number chemical potential, these other situations can 
be studied both analytically and numerically, and interesting new phases of matter have 
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been found there which should be investigated further. 
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TABLES 



11, 


T 

A ' I A 


T 


T 

^ A ' z A 


.25 


0617(2) 


0.1134(2) 


1962C4 


.28 


1120(5) 


1732(4) 


2546(5 


.29 


1407(5) 


1993(5) 

V / • J- \ > * J KJ \ KJ I 


2782(4 


.30 


1772(9) 


2269(6) 

v / • — — kj \ J \ v / j 


3016(5 


.31 


2120(6) 


2544(6) 


3240(3 


.32 


2415(11) 


2783(7) 

KJ * £J 1 KJ ' \ 1 / 


3470(6 


.33 


0.2711(7) 


0.3045(7) 


0.3687(5 


.34 


0.2951(11) 


0.3253(9) 


0.3893(6 


.35 


0.3200(9) 


0.3499(7) 


0.4082(7 


.36 


0.3437(12) 


0.3651(9) 


0.4287(6 


.38 


0.3815(13) 


0.4044(9) 


0.4627(7 



TABLE I. Diquark Condensates on a 16 4 lattice at A = 0.0025, 0.005, and 0.010 in the second, 
third, and fourth column, respectively. The first column lists the \i values. 
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< X T T2X > 



< X T T2X > 



.25 


0.0100(2) 


0.0031(2) 


.28 


0.0508(5) 


0.0371(5) 


.29 


0.0821(5) 


0.0693(5) 


.30 


0.1275(9) 


0.1193(9) 


.31 


0.1710(6) 


0.1664(6) 


.32 


0.2050(11) 


0.2031(11) 


.33 


0.2380(7) 


0.2376(7) 


.34 


0.2650(11) 


0.2661(11) 


.35 


0.2900(9) 


0.2896(9) 


.36 


0.3220(12) 


0.3292(12) 


.38 


0.3590(13) 


0.3628(12) 



TABLE II. Diquark Condensates on a 16 4 lattice at A = 0.00. The first column lists the 
fj, values, the second records the linear extrapolations of the diquark condensates and the third 
records the quadratic extrapolations. 
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A* 


< XX > 


< XX > 


< XX > 


.25 


0.4123(5) 


0.4012(5) 


0.3887(5 


.28 


0.4025(8) 


0.3850(6) 


0.3697(6 


.29 


0.3934(7) 


0.3758(5) 


0.3616(4 


.30 


0.3816(9) 


0.3664(7) 


0.3529(5 


.31 


0.3686(6) 


0.3541(6) 


0.3444(5 


.32 


0.3523(11) 


0.3409(8) 


0.3335(6 


.33 


0.3381(6) 


0.3298(8) 


0.3234(5 


.34 


0.3242(11) 


0.3157(8) 


0.3125(6 


.35 


0.3091(10) 


0.3057(7) 


0.3026(7 


.36 


0.2987(12) 


0.2908(9) 


0.2920(7 


.38 


0.2733(12) 


0.2722(9) 


0.2711(7 



TABLE III. Quark-antiquark Condensates on a 16 4 lattice at A = 0.0025, 0.005, and 0.010 in 
the second, third, and fourth column, respectively. The first column records the \i values. 



A < X T T2X > < X T T~2X > < X T T2X > < X T T2X > 

.002 0.1239(6) 0.1315(4) 0.1678(8) 0.2194(14) 

.004 0.1678(4) 0.1732(3) 0.1890(5) 0.2034(7) 

.006 0.2031(4) 0.2073(2) 0.2190(4) 0.2308(5) 

.008 0.2303(4) 0.2338(2) 0.2440(4) 0.2566(4) 

.010 0.2526(3) 0.2579(2) 0.2662(3) 0.2753(4) 



TABLE IV. Diquark Condensates on a 16 4 lattice at four estimates of the critical fi, 0.2855, 
0.2870, 0.2920, and 0.2970 in the second, third, fourth, and fifth column, respectively. The first 
column records the A values. 
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A < XX > <XX> <XX> <XX> 

.002 0.4244(9) 0.4261(6) 0.4362(9) 0.452(2) 

.004 0.3882(6) 0.3875(4) 0.3836(5) 0.3780(9) 

.006 0.3730(4) 0.3726(3) 0.3668(4) 0.3619(5) 

.008 0.3622(4) 0.3611(3) 0.3566(4) 0.3511(4) 

.010 0.3525(4) 0.3528(2) 0.3464(3) 0.3415(4) 



TABLE V. Quark-antiquark Condensates on a 16 4 lattice at four estimates of the critical /i, 
0.2855, 0.2870, 0.2920, and 0.2970 in the second, third, fourth, and fifth column, respectively. The 
first column records the A values. 



Form of scaling function 


a + bx 


a + bx s 


a + bx + cx 2 / 3 


Number of parameters 


5 


6 


6 


He 


0.278(7) 


0.285(7) 


0.2948(7) 


fimag 


0.5(1) 


0.51(8) 


0.58(4) 


5 


1.1(2) 


1.6(4) 


2.28(3) 


X 2 /dof 


86 


74 


17 



TABLE VI. Determination of the critical chemical potential and two critical exponents using the static 
scaling hypothesis assuming different forms for the scaling function f(x) and using the data for diquark 
sources 0.0025 < A < 0.010, and 0.297 < [i < 0.340, i.e. 19 data points. The numbers of parameters and a 
measure of the fit quality are given for each fit. 
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